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A phenomenological analysis is presented of recent two-particle angular correlation data on rela- 
tive pseudorapidity (jj) and azimuth reported by the Compact Muon Solenoid (CMS) Collaboration 
for y/s = 7 TeV proton-proton collisions. The data are described with an empirical jet-like model 
developed for similar angular correlation measurements obtained from heavy ion collisions at the 
Relativistic Heavy Ion Collider (RHIC). The same-side (small relative azimuth), ^-extended cor- 
relation structure, referred to as the ridge, is compared with three phenomenological correlation 
structures suggested by theoretical analysis. These include additional angular correlations due to 
soft gluon radiation in 2 — > 3 partonic processes, a one-dimensional same-side correlation ridge on 
azimuth motivated for example by color-glass condensate models, and an azimuth quadrupole sim- 
ilar to that required to describe heavy ion angular correlations. The quadrupole model provides 
the best overall description of the CMS data, including the ridge, based on x 2 minimization in 
agreement with previous studies. Implications of these results with respect to possible mechanisms 
for producing the CMS same-side correlation ridge are discussed. 

PACS numbers: 12.38.Qk, 13.85.Hd, 25.75.Bh, 25.75.Gz 



I. INTRODUCTION 

The CMS Collaboration at the Large Hadron Collider 
(LHC) at CERN recently reported the observation of 
long-range two-particle angular correlations on relative 
pseudorapidity (rj a = rji — 772) for charged particle pairs 
at mid-rapidity with relative azimuth (0a = 4>i — $2) 
< tt/2 (near-side or same-side - SS) in high multiplicity 
proton-proton (p-p) collisions at y/s = 7 TeV [l|. The 
collaboration noted the visual similarity of this same-side 
"ridge" with two-particle angular correlation structures 
observed in heavy ion collisions at RHIC For the 

latter, two-dimensional (2D) angular correlation analy- 
sis reveals at least two structures which produce same- 
side, relative 77 elongated correlations: (1) a quadrupole, 
cos 20a j independent of r\ (within two units at mid- 
rapidity) whose amplitude is approximately proportional 
to the square of the eccentricity of the transverse dis- 
tribution of colliding nucleons from the impacting nu- 
clei @, and (2) a 2D peak centered at (?7a,0a) = (0,0) 
whose amplitude and width on relative pseudorapidity in- 
crease markedly with more central collisions (i.e. reduced 
impact parameter) The quadrupole correlation in 
heavy ion collisions is conventionally attributed to pres- 
sure driven, hydrodynamic flow Q due to rapid equilibra- 
tion of a partonic medium ||. Inconsistencies between 
this assumed scenario and other RHIC data have been 
found however Q; alternate mechanisms f or q uadrupole 
correlation production have been studied [lOL [Tl|. The 
origin of the 77-elongated 2D peak at RHIC is not un- 
derstood although many ideas have been proposed. The 
appearance of one or both of these ridge-like correlation 
structures in p-p collisions at the LHC would be quite 
interesting. 

A number of mechanisms have been proposed to ex- 
plain the 77-elongated 2D correlation peak observed in 



the heavy ion collision data at RHIC. Several of these 
authors have extended their models to include the CMS 
p-p results. Dumitru et al . II 211 proposed that flux tubes 
of color glass condensate [13j formed immediately after 
impact produce correlations among the hadron fragments 
which arc transported outwar d by an assumed radial 
flow. Shuryak [Tj] and Voloshin [15| invoked similar ideas 
except the initial stage correlations result from beam-jets, 
which are the fragmentation results of forward going nu- 
cleon remnants following semi-hard partonic scattering. 
Hwa et al. [la ] proposed that the same-side ridge is gen- 
erated by fast parton (from semi-hard scattering) - soft, 
"thermal" parton recombination into final-state hadrons. 
Wong [17] explained the same-side ridge by way of scat- 
tering (momentum kick) between fast, transverse scat- 
tered partons and soft partons of a medium. Levin and 
Rezaeian [l8[ argued that the long-range rapidity and 
azimuth correlations between charged hadron pairs pro- 
duced from two BFKL parton showers will contribute to 
the CMS ridge. Werner et al. [H| interpreted the CMS 
results in terms of energy density fluctuations in the ini- 
tial collision stage which manifest in the final-state via 
hydrodynamic expansion; the relevant correlation being 
a quadrupole. 

Within perturbative QCD Field |20j suggested that ini- 
tial and final state radiation effects in partonic 2 — > 3 and 
2 — > 4 processes may contribute to extending the n range 
of the jet correlations. Sjostrand [2l[ and Trainor [22[ in- 
dependently suggested that fragmenting color flux tubes, 
as in the LUND model of soft particle production p3| . 
stretching between transverse scattered colored partons 
and longitudinal colored nucleon remnants from semi- 
hard collisions could produce 77 elongation in the jet cor- 
relation peak. None of the preceding models have been 
shown to account for all of the correlation 24| and pt 
spectrum [2|| features in the heavy-ion data from RHIC 
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which are associated with the 77-elongated peak phe- 
nomenon. 

In a physical model-independent analysis Trainor and 
Kettler [26| showed that the systematic dependences of 
the quadrupole correlation amplitude on collision central- 
ity, event multiplicity, and collision energy (~ log y/s) 
determined from RHIC heavy- ion data [6|], when ex- 
trapolated to LHC energies for high multiplicity p-p 
collision events, predicted quadrupole amplitudes which 
were in reasonable agreement with the CMS near-side 
ridge. They suggested that the CMS ridge and the 
quadrupole correlations at RHIC have a common dynam- 
ical origin, but not necessarily hydrodynamic. Bozek (27j 
also showed that the transverse momentum (p t ) se- 
lected azimuth projections of the CMS 2D correlations at 
larger relative pseudorapidity [l[ were well described with 
cos 4>a (dipole) and cos 20a (quadrupole) terms only. 

In the present analysis the 7 TeV p-p angular corre- 
lation data from the CMS experiment were fitted in 2D 
angular space (t7A)<Aa) using a multi-component model 
similar to that used to describe angular correlations from 
Au-Au collisions at RHIC [3]. Four data sets were 
studied corresponding to combinations of event selection 
(minimum-bias or high multiplicity) and particle pt selec- 
tion {p t > 0.1 GeV/c or 1 < p t < 3 GeV/c). Additional 
components (soft gluon radiation, azimuth quadrupole 
and same-side Gaussian ridge on azimuth) were included 
in the fitting model in order to phenomenologically rep- 
resent the correlation structures expected for each class 
of theoretical models listed above. The purpose of this 
analysis is to determine which phenomenological descrip- 
tion^) of the p-p CMS correlation data is preferred based 
on fit quality in an effort to guide ongoing theoretical 
studies of the CMS ridge. 

The paper is organized as follows. The phenomenolog- 
ical fitting model is explained in the next section; fitting 
results are presented in Sec. lIIII The relation of this anal- 
ysis to previous studies of these data and the implications 
of the fitting phenomenology for theoretical models of rj- 
elongated angular correlations in high energy collisions 
are discussed in Sec. IIVI A summary and conclusions are 
given in Sec. |Vj 



taking the two particles in a pair from different but sim- 
ilar events. If events produce more than one jet within 
the acceptance all intra-jet charged particle pairs con- 
tribute to the correlation peak at small relative angles 
near (0,0). Inter-jet two-particle correlations are indistin- 
guishable from background for uncorrelated jets. Particle 
pairs from separate jets of a correlated, back-to-back dijet 
appear as a broad, opposite-side or away-side (4>a > tt/2, 
AS) ridge on azimuth centered at w. The same-side rj- 
elongated correlations may or may not be associated with 
jets; proximity of correlation structures on relative angle 
does not necessarily imply proximity in absolute (rj, 4>) 
space. 



B. Nominal fitting model components 

The primary elements of the present fitting function 
were developed previously by the STAR Collaboration in 
order to describe 2D angular correlation data for charged 
particle production from minimum-bias p-p collisions at 
y/s = 0.2 TeV [28|. Those studies used pt cuts and like- 
sign and unlike-sign pair combinations to reveal simple 
geometrical shapes in the correlations. Those structures 
include: (1) a jet-like, 2D Gaussian peak at small relative 
angles with an accompanying away-side, r?A-independent 
ridge corresponding to approximately back-to-back di- 
jets (e.g. pythia predictions in |lj) or other p t con- 
serving correlation effects, (2) a narrower 2D exponential 
peak at zero opening angle representing quantum correla- 
tions [29[ among identical charged particles, and (3) a ID 
Gaussian ridge on 77A which is the 2D manifestation of 
charge-ordering first observed on relative rj in collisions at 
the ISR [30| . The latter is thought to represent longitu- 
dinal (beam direction) fragmentation as described for ex- 
ample in the LUND model [23[ . With the addition of an 
^-independent quadrupole this phenomenological model 
quantitatively describes the 2D angular correlation data 
from minimum- bias Au-Au collisions at STAR [li]. 

C. Color coherence in 2 — » 3 processes 



II. ANGULAR CORRELATIONS AND FITTING 
MODEL 

A. Minimum-bias angular correlations 

The two-particle angular correlation distributions from 
CMS [1] and the STAR experiment at RHIC 2 4] repre- 
sent inclusive sums of all charged particle pairs within the 
tracking acceptances of the respective detectors. There 
is no special leading particle or "trigger" particle in this 
type of analysis and in this sense the CMS angular cor- 
relations may be referred to as minimum-bias. Corre- 
lations are obtained from normalized, binned ratios of 
real to mixed-event pairs, the latter being constructed by 



In addition to the above model components important 
terms could arise from color coherence effects [3l|, |32| . 
Soft gluon radiation amplitudes from color connected 
partons in hard scattering 2 — > 2 processes interfere pro- 
ducing increased yields in the hard scattering plane de- 
fined by the incident and scattered partons. Assuming 
hadronization does not obliterate the soft gluon distri- 
bution the resulting final state hadrons will increase the 
particle yield near the accompanying jet with a larger 
increase occurring along the beam direction, or on rela- 
tive pseudorapidity than on relative azimuth. This next- 
to-leading-order process may therefore contribute to the 
CMS ridge. For dijets within the acceptance soft gluon 
radiation will also contribute to the away-side correlation 
ridge on azimuth. 
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Color coherence effects can be approximated within 
the soft gluon momemtum limit (k — > 0) ;31| where the 
squared amplitude for the 2 — > 3 process is factored into 
the product of the probability density for the 2 — > 2 hard 
scattering times the soft gluon radiation probability dis- 
tributions for each color connected parton pair. The to- 
tal yield is obtained by summing over all color connected 
parton pairs. The dominant 2 — > 3 reaction considered 
here for midrapidity hadron production is gluon-gluon 
scattering given by [3l| 

g{pi)+g{p2) -> g{p3) + g(p4)+g{k), (l) 

where the soft gluon yield is approximately 

4 

P g « C A Y, Wa, (2) 

i<j=l 

the summation includes all parton pairs, gluon color fac- 
tor Ca — N c (number of colors) , and 



In Eq. ([3]) g is the strong coupling constant and pi -pj etc. 
denote scalar products of 4-momentum vectors. Quan- 
tities Wij are singular for massless partons in the limits 
k — > 0, 9ik —> 0, or 9jk — > where 6 are the angles be- 
tween the parton and radiated soft gluon. To regulate 
the singularities small cut-off masses (5 MeV/c 2 ) were 
assumed for each parton and a Gaussian angular cut- 
off, Ft = [1 - exp(— # 2 fc /2cr 2 )] where a g = 0.2 radians, 
was applied to each amplitude and Wij was weighted 
with FiFjWij. The latter cut-off also respects angular 
orde ring requirements for the subsequent jet fragmenta- 
tion [31| • The soft gluon radiation yield in the soft limit 



was approximated by 

where — 1 — /3ij3j cos %, (3 = p/E (velocity), gluon 
1 defines the beam direction, Ex k is the transverse en- 
ergy of the radiated gluon, and A g is a phenomenologi- 
cal amplitude which will be adjusted by the fitting pro- 
gram. The final-state hadron distribution resulting from 
the soft radiated gluons is assumed to be given by Eq. ^ 
accordin g to the local parton-hadron duality hypothesis 
(LPHD) [Hj]. 

The soft gluon radiation adds to the final-state hadron 
distribution associated with minimum-bias jet produc- 
tion, given by Pj c t,3 + -Pjct.4 + Pg- The first two yields 
correspond to hadron distributions from the fragmenta- 
tion of hard scattered partons 3 and 4 in the reaction in 
Eq. (P). Each is represented by a 2D Gaussian distribu- 
tion 

P ]ct { v ,4>) = Ae- (j) -^ t)2/2(T 'e^ w "^ t)2/2CT *, (5) 
where 77, </> are defined within the detector acceptance. 
The single particle Gaussian widths are directly propor- 
tional to the widths of the jet-like angular correlation on 
77A and 4>a discussed below. It is expected that P g <C Pj e t 
and the additional contribution to the 2D angular corre- 
lations is dominated by the cross term P g (Pj e t,3 + -Pjet, 4) 
projected onto relative azimuth and pseudorapidity and 
averaged over all jet angular positions within the accep- 
tance. The 2D angular autocorrelation in the soft gluon 
radiation limit for 2 — > 3 processes was approximated by 
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where Arj is the (symmetric) pseudorapidity acceptance, 
0(t7a) = |Ary-|ry A ||, Vt, = m+m and S = 0i+0 2 . Jets 
3 and 4 were for simplicity assumed to be uniformly dis- 
tributed within the 2ttAti acceptance and back-to-back 
in the p-p collision cm. system (CMS lab frame), where 

</>jct,4 = 0jct,3 ± 7T and 77jet,4 = ~?7jet,3- Pg(v^) was ob- 
tained from Eq. (j4]) for back-to-back dijets at angular po- 
sitions (j7j e t,3, 0jct,3) and (?7jet,4, <#jet,4) with the soft gluon 
(resulting final-state hadron) at (rj,(p). 

The resulting angular correlations from Eq. ([6]) for the 
high multiplicity, 1 < p t < 3 GeV/c data set are shown 
in Fig. [T] for different values of cut-off parameters and for 



jet widths er j; = 0.20 and — 0.22. The parton and soft 
gluon cut-off masses and collinear cut-off a g assumed in 
panels (a) - (d) were 5 MeV and 0.2 radians, 1 MeV and 
0.2 radians, 5 MeV and 0.05 radians, and 5 MeV and 0.35 
radians, respectively. In each panel the amplitude was 
arbitrarily adjusted to unit height difference between the 
positions on (t]a,4>a) at (0,0) and (0,7r/2). The away- 
side ridge corresponds to correlations between soft-gluons 
emitted near jet-3 and the final-state hadrons from jet-4 
and vice versa. The increased amplitude of the away-side 
ridge at larger |?7a| results from the increased radiative 
yields when the angles between incident and outgoing 
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(a) (b) 




FIG. 1: (Color online) Perspective views of two-dimensional 
charge- independent angular correlations from Eq. (O, as- 
suming arbitrary normalization, for soft-gluon radiation from 
back-to-back Gaussian jets estimated in the soft limit (k — > 0) 
on relative azimuth and pseudorapidity (??a,<^a) f° r the high 
multiplicity, 1 < p t < 3 GeV/c data from CMS. From upper- 
left, to lower-right panels (a)-(d) show angular correlations 
assuming nominal cut-off masses 5 MeV and angle 0.2 rad; 
1 MeV, 0.2 rad; 5 MeV, 0.05 rad; and 5 MeV, 0.35 rad, re- 
spectively. The single particle jet distribution widths on rj 
and (ft were 0.20 and 0.22 respectively. 



particle yields from relativistic heavy ion collisions is of- 
ten cited as evidence for hydrodynamic evolution, other 
theoretical studies [l(| EH show that a quadrupole cor- 
relation may result directly from QCD. 

Models assuming outward, collective flow with initial- 
state energy fluctuations {e.g. CGC, beam-jets) or re- 
combination suggest including a same-side ID ridge on 
azimuth with no accompanying away-side structure other 
than that already present in the model. Momentum con- 
servation induced correlations associated with this class 
of models would manifest as an away-side ridge on rjA and 
are therefore already included in the fitting function. 



E. Specific modifications for the CMS data 

Other modifications of the fitting model in [3j were 
required for the CMS data. At RHIC energies the away- 
side correlation ridge (non-quadrupole) for lower p t par- 
ticles is well described by a dipole, cos 0a, which is un- 
derstood to represent the limit of a periodic away-side 
azimuth ridge (Gaussian) as its width increases . Pre- 
liminary studies showed that the 7 TeV CMS correlation 
data were better described with a periodic Gaussian. The 
measured correlations also display significant away-side 
?7A-dependence. Analysis of Au-Au quadrupole correla- 
tions with large r\ acceptance indicates that reduction in 
amplitude with increasing |?7a| is possible (36| . Conse- 
quently, the away-side ID Gaussian, the same-side ID 
Gaussian, and the quadrupole components all include 
symmetric r^A-depcndent modulations. 



partons decrease. The fall-off at the largest \t]a\ results 
from both the collinear cut-off a g and the sin 2 9ik factor 
in Eq.©. 

The angular correlations are weakly sensitive to the 
assumed cut-off mass (panels a and b), but for masses 
of order the pion mass the correlations become similar 
to that in panel (d). The angular correlations are how- 
ever strongly affected by the collinear cut-off. For the 
fits to data a g = 0.2 was selected (panel a) in considera- 
tion of angular ordering while minimizing the suppression 
at larger |77a|- For the fits to the four 7 TeV data sets 
cut-off parameters 5 MeV and 0.2 radians were always 
used. The jet single particle Gaussian widths on rj and (ft 
were (0.35,0.52), (0.22,0.26), (0.30,0.30) and (0.20,0.22) 
for the minimum-bias low and higher p t data, and the 
high multiplicity, low pt and higher p t data sets, respec- 
tively. 



D. Ridge components 



Recent theoretical predictions [19|,|34| and phenomeno- 
logical analysis of the 7 TeV data from CMS [H, H3| 
warrant including the quadrupole component. Although 
the presence of the quadrupole correlation in the charged 



F. Fitting model 

The complete fitting function used here is given by [37] ] , 
F(r]A,(ftA) = A + A Q f{n A ) cos(20 A ) 
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k—even—int 



where /(?7a) = 1 + Sr] A , g(r]A) = 1 
C cos(27t?7a/A?7), and the summations for the peri- 
odic away-side and same-side ID Gaussian ridges (A4 
and A§ terms) were accurately truncated with sums 
(-3,-1 •••5) and (-4, -2 •••4), respectively. The 
quadrupole {Aq term) and the same-side periodic ID 
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Gaussian (A$ term) were applied alternately, not simul- 
taneously in the fitting. Terms periodic in </>a for the Ai 
and A2 components [i.e. dependent on 0a — (±27r)] were 
included but were negligible for the azimuth widths used 
here. The constant offsets for the AS and SS Gaussian 
ridge components were included in order that coefficients 
Ai and A§ only (approximately) affect the amplitudes of 
the </>A-dependent oscillations and not the overall offset. 

Function F(rj was fitted to measured quantity 
-R^a^a) from CMS [l|. The fitting parameters in- 
clude: A (normalization offset), Aq or A 5 and ass with 
5, same-side jet parameters A\, cr VA , o<p A , a and /3, 2D 
exponential parameters A%, and w<f,, ID Gaussian pa- 
rameters A3, (To and 7 (only required for the minimum- 
bias, pt > 0.1 GeV/c data), and away-side ID Gaussian 
ridge parameters A4, a as, £ and £. Due to the similar 
magnitudes of the width parameters for the A\ and A2 
terms the fitting model may not provide an accurate sep- 
aration between jet fragment correlations and quantum 
correlations in each case. 

The soft gluon radiation component was computed for 
each case using the Gaussian width parameters for the 
jets obtained by fitting the data without the soft radi- 
ation component and with the cut-off parameters dis- 
cussed in the above subsection. The single particle jet 
width parameters in Eq. ((SJ are related to the SS 2D 
correlation peak Gaussian widths by = <t, ;a /y/2 and 
00 = <T0 A /v / 2 owing to the definition of the difference 
variables 77A and 0a- The function computed in Eq. ^ 
was arbitrarily adjusted in overall magnitude by the fit- 
ting program in order to obtain the best fit. 



G. Fitting errors 

Best fits were selected by x 2 P er degree of freedom 
(dof), residuals, the 2D x 2 distributions, and parameter 
fit errors. The x 2 /dof values were fairly large indicating 
that residual, systematic effects are statistically signif- 
icant. Such effects may include experimental artifacts 
and/or physical processes which are not accurately de- 
scribed by the fitting model. The fit parameter error 
covariance matrix was estimated by the following matrix 
inversion, 



S 2 P 



3* 



E 



1 dF t 8F t 
e\ dPj dP k 



(8) 



where the summation includes all unique (t?a,0a) bins, 
e, is the statistical error, and dFi/dPj are partial deriva- 
tives of the model function with respect to fit parameter 
Pj at the x 2 minimum. This estimate assumes that near 
the x 2 minimum the fitting function dependence on pa- 
rameter Pj can be accurately expanded to the leading 
order (linear) terms in a Taylor series. Parameter fit er- 
rors, V 5 2 P ii , are listed along with the best fit values in 
Table |U The CMS collaboration quotes an overall 15% 
systematic uncertainty in the data near the same-side 



ridge resulting in an additional 15% systematic uncer- 
tainty in the quadrupole and same-side ridge amplitudes 
Aq and A$. 



III. RESULTS 

The fitting results for the high event multiplicity data 
with 1 < p t < 3 GeV/c from the CMS experiment are 
discussed first. Initially the data were fit using the Ao, 
Ai, A 2 and A4 terms in Eq. ([7]) where ID Gaussian 
component A3 was not required. The same-side peak 
structure was accurately described by the combination 
of approximate 2D Gaussian and 2D exponential terms 
where the fitted exponents of the former term equaled 2 
(Gaussian) within 2% (but note the caveat in Sec. Ill Fl) . 
The away-side ridge was well described by the periodic 
Gaussian (A4 term) with moderate tja modulation. Al- 
though the SS peak and AS ridge were well described 
the x 2 was rather large due to the inability of the model 
to account for the ^-extended same-side ridge. Includ- 
ing the soft gluon radiation correlations with adjustable, 
non-negative amplitude failed to improve the fit qual- 
ity in this instance. x 2 /dof was improved significantly 
however (by 12 units) when either the quadrupole or SS 
periodic Gaussian ridge (A5 term) was included. Fits 
including the quadrupole obtained the lowest x 2 ( see Ta- 
ble [3J. The data, the best fit (with quadrupole), and the 
residuals (model - data) are shown in Fig. [5] in the right- 
most column of panels. The soft gluon radiation corre- 
lations afforded very little improvement in fit quality for 
the SS Gaussian ridge model, but significantly improved 
the x 2 /dof from 4.18 to 4.04 when the quadrupole was 
included. 

Similar results were found for the high multiplicity, 
Pt > 0.1 GeV/c data except that the x 2 improvement 
obtained with either the additional quadrupole or SS 
Gaussian was less dramatic (about 4 units for % 2 /dof) 
than was obtained for the higher p t data. The SS 2D 
peak function (A±) relaxed to a non-Gaussian geom- 
etry with both pseudorapidity and azimuth exponents 
less than 1 (see Table [I]). Statistically significant evi- 
dence for both the quadrupole and SS Gaussian ridge 
was found. The best fit was obtained assuming the addi- 
tional quadrupole component. The soft gluon radiation 
correlations provided no significant improvement in the 
fit quality. The data, best fit (quadrupole model), and 
residuals are shown in Fig.[5]in the third column of panels 
from the left. 

The quadrupole and same-side Gaussian ridge can have 
very similar shapes and amplitudes for the same-side az- 
imuth range. The two are of course very different on the 
away-side. Evidently the % 2 /dof reduction afforded by 
the quadrupole for both high event multiplicity data sets 
is due, at least in part, to improvements in the fit qual- 
ity for the away-side data. Visual evidence of an away- 
side quadrupole contribution to the CMS high multiplic- 
ity data is not apparent, unlike the angular correlation 
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data for mid-central relativistic heavy ion collisions [3j,l6|. 
Quantitative analysis of the type presented here and in 
Refs. [H, H3] is required to differentiate between these 
two descriptions of the SS and AS ridges. Both the 
quadrupole and SS Gaussian amplitudes increase for the 
higher p t selected particle cuts. A modest, negative T)\ 
modulation was obtained for both the quadrupole and 
SS Gaussian which had only minor effects. 

The magnitudes of the residuals for the high event mul- 
tiplicity data are less than 0.1 except near the origin and 
the ?7a acceptance edge and are about 1% of the mag- 
nitudes of the principal correlation structures. Although 
the residuals increase at large |?7a| the increasing statis- 
tical errors at large | 77 a | prevent these bins from making 
significant \ 2 contributions. The band structures in the 
residuals on 77a and </>a are not understood and could 
be of experimental and/or dynamical origin. In any case 
the residuals are small relative to the principal correla- 
tion structures including the quadrupole and SS Gaussian 
ridge. 

Both sets of minimum-bias correlation data were well 
described with just the Ao, A\, A%, A3 (only for the 
Pt > 0.1 GeV/c data) and A4 components, x 2 f° r the 
Pt > 0.1 GeV/c data was dominated by the single (0,0) 
angular bin which was subsequently omitted in the \ 2 
minimization procedure. No evidence for a non-negative 
quadrupole or SS Gaussian ridge structure was found. 
However, including the soft gluon radiation correlations 
improved the x 2 /dof by about 1 unit. The data, best fits 
(no quadrupole or SS Gaussian ridge), and residuals are 
shown in the first two columns of panels in Fig. [21 The 
magnitude of the residuals is less than 0.03 except near 
the origin and the 77 a acceptance edge and is about 1% 
of the magnitudes of the principal correlation structures. 

For each case the same-side peak structure was well de- 
scribed with a combination of approximate 2D Gaussian 
and 2D exponential functions. The exponent fit param- 
eters a and /3 were generally close to unity correspond- 
ing to a 2D Gaussian-like peak shape. For each data 
set the five parameters of the 2D Gaussian-like compo- 
nent were not significantly affected by including or not 
including the quadrupole and SS Gaussian ridge. Trainor 
and Kettler [26| showed that the minimum-bias jet cor- 
relation structure in the CMS data is similar to that at 
RHIC energies when scaled by log(i/s). The total num- 
ber of jet-like correlated pairs per final-state particle, es- 
timated by the volume of the same-side 2D Gaussian-like 
peak (see Table HI but note the caveat in Sec. Ill F|) . in- 
creases for high multiplicity events because those events 
have increased probability of minimum-bias jet produc- 
tion within the acceptance. However, this volume de- 
creases for the higher p t (1 < pt < 3 GeV/c) particle 
selection cut. Studies of 2D charged particle (pti,pt2) 
correlations from y/s = 0.2 TeV minimum-bias NSD p-p 
collisions at RHIC [H] suggest that the p t £ [1, 3] GeV/c 
cut-range may exclude a large fraction of the correlated 
pairs produced by minimum-bias jets. This may account 
for some of the reduction. The widths of the same-side 



jet-like peak decrease with higher p t selection as expected 
for jet fragmentation with uniform transverse momentum 
distributions relative to the jet thrust axis [38| . Further 
analysis of the minimum-bias jet-like correlations in p-p 
collisions at the LHC will likely require different pt selec- 
tion cuts which are optimized for that purpose. 

In general the away-side azimuth ridge structure was 
better described with an r/A-dependent, periodic Gaus- 
sian than with an ?7A-dependent dipole. A periodic se- 
ries of ID Gaussians, as in Eq. ([7]), approaches an az- 
imuth dipole plus offset for sufficiently large width a as 
{e.g. higher-order multipole contributions arc < 1% for 
a as > 1-75 Hl|). The fitted widths here are less than 
1.2 indicating that the dipole limit is not reached. The 
amplitude A4 approximately follows the volume trends 
of the SS 2D peak as expected if this correlation struc- 
ture is dominated by inter-jet charged particle fragment 
pairs for dijets within the tracking acceptance. The ma- 
jority of the ?7a dependence (fall-off) of the away-side 
ridge (see Fig. [5J is described by the negative e?? A mod- 
ulation. Weaker dependences were approximated by the 
Ccos(27T77a/A77) modulation. Model descriptions of the 
AS 7/a dependence did not significantly affect the fitted 
parameters of the SS 2D peak or the SS ridge structure. 
The away-side Gaussian decreased in amplitude and in- 
creased in width when the quadrupole component was 
included in the fitting. 



IV. DISCUSSION 

The CMS collaboration noted the similarity between 
the same-side 77A elongated correlations and that re- 
ported for Au-Au collisions at RHIC by the STAR @- 
|4j and PHOBOS [a] collaborations for total energy per 
nucleon-nucleon (N-N) pair ^snn — 0.2 TeV. What was 
observed in minimum-bias collisions at RHIC using all 
charged particles with \r]\ < 1.0 and pt > 0.15 GeV/c 
was a 2D approximately Gaussian peak centered at 
(?7A, 4>a) — (0, 0) whose amplitude and 77 width increase 
substantially and monotonically with decreasing impact 
parameter between the colliding ions (increasing central- 
ity) . In addition analysis of these data provided evidence 
for an rj^- independent quadrupole which smoothly varies 
with collision centrality and is approximately propor- 
tional to the square of the spatial eccentricity of the over- 
lapping, colliding ions. In similar correlation studies us- 
ing higher momentum "trigger" particles with lower mo- 
mentum associated particles the same-side peak evolved 
to a non-Gaussian structure [4|. It is plausible that jet 
fragmentation could be modified in the dense environ- 
ment produced in relativistic heavy ion collisions [39] . 

Field [2(| pointed out that initial and final state radia- 
tive processes (2 — > 3, 2 — > 4) accompanying transverse 
partonic scattering preferentially produce hadronic frag- 
ments in the hard-scattering plane defined by the beam 
and jet axes. Single gluon radiation in 2 — ¥ 2 hard pro- 
cesses was calculated in the soft gluon limit in [3l|, [32j . 
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MB p > 0.1 GeV/c MB p ( s [1 ,3] GeV/c N> 110, p>0.l GeV/c N> 110, p t e [1 ,3] GeV/c 






FIG. 2: (Color online) Perspective views of two-dimensional charge-independent correlations for p-p minimum-bias and high 
event multiplicity (N > 110) collisions at ^/s — 7 TeV from the CMS collaboration [l[. Upper, middle and lower rows of panels 
display data, model fits, and residuals (model - data), respectively. The columns of panels from left to right correspond to 
minimum-bias events with charged particle selection pt > 0.1 GeV/c and 1 < p t < 3 GeV/c, and high multiplicity events with 
pt > 0.1 GeV/c and 1 < pt < 3 GeV/c, respectively. The fitting models selected for display are described in the text. All 
residuals are shown using the same, expanded scale. 



Using the gg — > ggg distribution from (3l| additional an- 
gular correlations for the p-p collision system were cal- 
culated and included in the present fitting model. Al- 
though the additional correlations made statistically sig- 
nificant improvements in the fit quality for the minimum- 
bias collision data, they could not explain the same-side 
7?A-extended ridge in the high event multiplicity data. 



Sjostrand [2l| and Trainor [22| suggested that color- 
flux tube connections between transverse (scattered) and 
longitudinal (beam direction) partons could, in analogy 
to the LUND color-string fragmentation model [23| for 
soft particle production, produce an ?/a extended, ex- 
ponentially decreasing enhancement in the jet-like corre- 
lations. The present analysis does not find evidence of 
increasing rj width (i.e. increased value of a v& ) or evo- 
lution toward an exponential shape [exponent parameter 
a in Eq. ([?])] in conjunction with the appearance of the 
same-side ridge, in agreement with previous findings by 
Trainor and Kettler [26[. Nevertheless, this suggested jet 



fragmentation mechanism should be further considered 
with respect to the ?yA-elongated correlations observed 
in A- A collisions at RHIC and the LHC 0. 

Two classes of theoretical models purport to explain 
the same-side ^-elongated ridge correlations observed 
in the RHIC heavy ion data. The first assumes event-by- 
event energy and/or momentum density fluctuations in 
the initial collision stage immediately after impact which 
propagate outward via pressure driven radial expansion. 
The initial stage fluctuations are described as color-glass 
condensate (CGC) [l3| flux tubes (glasma) 01, localized 
"hot spots" in a quark-gluon plasma (T^ . l4lll42l. or beam- 
jet remnants from hard scattering processes [l4fl5|. The 
second class of models assumes hard, transverse parton 
— soft, longitudinal parton interactions via recombina- 
tion [16[ or strong rescattering [l7( to induce the tja- 
elongated correlations. 

So far, most proponents have only attempted to ad- 
dress a subset of the experimental features associated 
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TABLE I: Best fit model parameters with fitting statistical errors (in parentheses) for p-p minimum-bias and high event 
multiplicity (N > 110) collisions at </s = 7 TeV. 

Parameter MB p t > 0.1 GeV/c MB p t G [1,3] GcV/c N > UO~p t > 0.1 GeV/c N > 110, G [1,3] GcV/c 





Jct+ASG 


Jet+ASG 


Quad 


SSG 


Quad 


SSG 




1.915(19) 


4.40(3) 


11.63(15) 


11.90(17) 


16.52(15) 


16.17(14) 


°"<?A 


0.501(3) 


0.313(1) 


0.426(4) 


0.422(4) 


0.285(1) 


0.285(1) 


a 


1.021(5) 


1.032(4) 


0.938(8) 


0.916(7) 


0.979(5) 


0.990(5) 


°>A 


0.737(7) 


0.370(1) 


0.422(7) 


0.475(6) 


0.310(1) 


0.318(2) 


/8 


0.883(9) 


1.035(4) 


0.740(6) 


0.794(7) 


1.003(6) 


1.012(7) 


A 2 


8.07(3) 


6.88(3) 


35.0(2) 


34.8(2) 


14.3(2) 


14.68(15) 


w,, 


0.753(2) 


0.462(1) 


0.446(2) 


0.446(3) 


0.418(3) 


0.421(3) 




0.807(5) 


0.429(1) 


0.458(3) 


0.435(2) 


0.377(2) 


0.363(2) 


A 3 


1.532(24) 


— 


— 


— 


— 


— 


(To 


0.826(3) 


_ 


_ 


— 


— 


_ 


7 


0.931(3) 


— 


— 


— 


— 


— 


A 4 


2.125(12) 


1.446(1) 


5.62(2) 


6.23(3) 


3.08(2) 


3.80(2) 




1.196(2) 


0.956(1) 


0.947(8) 


0.862(4) 


1.156(31) 


0.841(5) 


e 


-0.021(1) 


-0.0326(1) 


-0.0231(3) 


-0.0242(3) 


-0.020(1) 


-0.0208(3) 


c 


-0.055(3) 


0.0080(6) 


-0.040(1) 


-0.035(1) 


-0.046(2) 


-0.042(1) 


Aq 






0.486(23) a 




0.657(24)° 












0.556(40) a 




0.846(27) a 










0.550(28) 




0.583(15) 


5 






-0.056(2) 


-0.017(6) 


-0.032(1) 


-0.021(3) 


A 


-1.344(5) 


-0.124(1) 


-1.038(18) 


-0.571(17) 


-0.505(19) 


0.164(11) 


X 2 /dof 


14.3 


16.8 


9.78 


10.9 


4.18 


4.57 


Vol(2DG) 


4.64 


3.13 


15.5 


17.2 


9.23 


9.20 




2.54 


1.38 


5.32 


5.37 


3.56 


3.20 


«2 






0.046 




0.053 





"Additional 15% systematic uncertainty should be included due 
to quoted systematic error in the ridge signal pj . 



with the same-side angular correlations, e.g. either the 
width or amplitude dependence on collision centrality. 
Accompanying correlation structures including the away- 
side ridge, charge-dependence (like - unlike charged pair 
difference), and the 2D two- particle (pn,Pt2) correla- 
tion |24[ have not yet been addressed by these theo- 
retical models. A recent analysis [43[ of [ptiiPti] cor- 
relation predictions for CGC flux tube induced correla- 
tions [44| shows that such initial stage models, even in- 
cluding strong radial flow, are irrelevant for describing 
the same-side r?A-extended correlations. Nevertheless, 
these models motivated the periodic, same-side Gaussian 
ridge component adopted in the present phenomenolog- 
ical study. Fits to the minimum-bias data provide no 
significant evidence for this component while fits to the 
high event multiplicity, N > 110 data do. 

The quadrupole angular correlation continues to be ex- 
tensively studied in relativistic heavy ion collisions Q 
and was recently reported for Pb-Pb collisions at the 
LHC by the ALICE collaboration (45[. The most com- 
pelling explanation for the quadrupole correlation is that 
it is a consequence of anisotropic particle emission of the 
form [1 + 2v2 cos(2(0 — ip))] , where ip is a random, event- 
wise angle for maximum azimuth particle density. Aver- 
aging over an event ensemble produces two-particle cor- 
relation component, 1v\ cos(20a)- The anisotropic emis- 
sion presumably reflects the initial geometry (i.e. ec- 



centricity e) of the transverse overlap region between 
the colliding constituents of the two nuclei. Measure- 
ments show that observed V2 is approximately propor- 
tional to initial eccentricity estimated using a Glauber 
model [(| . Naive superposition of independent N-N colli- 
sions, even if the spatial distribution of N-N collision po- 
sitions is non-isotropic, produces on average an isotropic 
final-state particle distribution on azimuth resulting in 
V2 = 0. What causes the azimuth modulation (v?) is 
a mystery although theories are abundant (e.g. Q but 
note discussion in Q). 



Trainor and Kettler [26[ and Bozek [27J both showed 
that the p-p high multiplicity 7 TeV CMS correlation 
data could be well described with a quadrupole compo- 
nent. Trainor and Kettler described the 7 TeV CMS 
data with a similar, 2D fitting model as in Eq. ([7]) but 
without the 2D exponential, same-side Gaussian ridge, 
and 77A modulations assumed here. They showed that 
the quadrupole amplitude inferred for the 7 TeV high 
multiplicity p-p correlation data is consistent with the 
simple, linear trends on collision energy [proportional 
to log(y/ snn)], eccentricity and event multiplicity re- 
ported for Au-Au minimum-bias collisions at RHIC [6]. 
Bozek [27} fitted the ID azimuth projections for |t7a| from 
2.0 to 4.8 for sixteen combinations of event multiplic- 
ity and p t interval [l[ using cos(<^a) and cos(2</>a) terms 
only. Good fits were obtained although contributions of 
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the same-side jet peak for \t)a\ > 2 occur and the azimuth 
projection averages over the strong 7/a dependence of the 
away-side ridge beginning at \t]a\ > 2. The quadrupole 
amplitude Aq in Eq. J?§ for the CMS correlation quan- 
tity R (see Eq. (1) in 1]) is related to azimuth anisotropy 
parameter v 2 defined above by, 



required. Therefore, if the p-p CMS correlation data 
do in fact display a quadrupole correlation as concluded 



offline 
irk 



}-l V, 



(9) 



Inferred values of v 2 from the present analysis (see Ta- 
ble HI and Eq. © are v 2 = 0.046 and 0.053 for the 
N > 110, p t > 0.1 GeV/c and 1< p t < 3 GeV/c data, re- 
spectively, in fair agreement with ~ 0.04 obtained in both 
Trainor and Kettler [2f| and Bozek [13|. The best fits 
obtained in this analysis were those with the quadrupole 
component as discussed in Sec. Mil 

The CMS experiment also reported strong variation in 
the large \t)a\ (2-0 < \t]a\ < 4.8) same-side ridge am- 
plitude for high multiplicity events (N > 90) in four 
Pt ranges from 0.1 to 4 GeV/c (see Fig. 9 in [l|). The 
larger values occurred in the pt € [1, 2] and [2, 3] GeV/c 
bins. A similar increasing then decreasing trend with p t 
was reported by the STAR experiment Q for the Au- 
Au quadrupole correlation in more central Au-Au colli- 
sions. CMS also found that the same-side ridge correla- 
tions for like-sign and unlike-sign charged particle pairs 
were the same within uncertainties as shown in Fig. 10 of 
[l| . Quadrupole correlations in Au-Au collisions at RHIC 
are also charge- independent within errors 0, H, @| . Non- 
zero charge-dependent angular correlations in Au-Au col- 
lisions at RHIC extend out to only 2 units on relative 
r] 0. The CMS results in Figs. 9 and 10 of Ref. Q are 
qualitatively consistent with quadrupole phenomenology 
for Au-Au collisions at RHIC energies. 

The conventional assumption in the heavy ion physics 
literature is that strong interactions among the (par- 
tonic) constituents lead to rapid thermal equilibration 
within about 1 fm/c after impact Q. This results in 
pressure gradients which drive collective expansion as 
described by particle transport models 47] or hydrody- 
namic evolution 0] ■ Hadronization models are invoked in 
the final stages of expansion which enable predictions to 
be compared with experiment. Event-wise fluctuations 
in energy and momentum density in the initial stage of 
the collision are expected to propagate to the final state 
producing correlations. Critical to this interpretation is 
the requirement of very rapid thermal equilibration; de- 
lays beyond a few fm/c result in insufficient hydrody- 
namic pressure to account for observed v 2 values. It has 
been shown [48| that rapid thermal equilibration cannot 
occur via pQCD processes or by non-perturbative instan- 
ton mechanisms in the brief amount of time available (of 
order 10 fm/c) during the heavy ion collision process. 
That thermal equilibration might occur during a p-p col- 
lision is even more problematic due to the reduced spatial 
and temporal scales involved. Kovchegov [48] concluded 
that if rapid thermal equilibration does occur some as 
yet unidentified, non-perturbative mechanism would be 



in [26J, [27J and indicated in the present results via \ se_ 
lcction, a production mechanism for the quadrupole cor- 
relation other than rapid thermal equilibration is implied. 

Three alternate mechanisms to the hydrodynamic sce- 
nario for generating the azimuth modulation in heavy ion 
collisions have been proposed. Boreskov et al. [lfj as- 
sumed a non-local N-N inclusive collision vertex model 
with final-state emission oriented relative to the mo- 
mentum transfer. Convolution of this vertex function 
with the transverse nucleus-nucleus overlap distribu- 
tion, which includes important anisotropic density gra- 
dients, produced anisotropic final-state particle distri- 
butions with tio values of reasonable magnitude. Kope- 
liovich et al. [ljj used a color-dipole model and predicted 
v 2 for p-p collisions at 200 GeV of order a few percent, 
which is a reasonable magnitude given the present results 
and those in [26], 27y and the empirical energy dependence 
of the quadrupole correlation reported by Kettler @. 
Trainor suggested that interacting color currents in 
the nucleus-nucleus overlap region could emit gluon mul- 
tipole radiation in analogy to classical electrodynamics. 



V. SUMMARY AND CONCLUSIONS 

Theoretically motivated phenomenological fitting 
models were applied to the two-particle 2D angular cor- 
relation data for p-p collisions at yfs = 7 TeV reported 
by the CMS Collaboration. The dominant features of 
the data were well described by jet-like correlation struc- 
tures projected onto relative 77, <j) subspace, consisting of 
a 2D Gaussian-like peak at the origin with accompanying 
away-side ridge at a = representing pair correlations 
between dijets within the acceptance. For the high multi- 
plicity events with N > 110 a new, extended correlation 
ridge along r]A was observed at small relative azimuth. 
The same-side jet-like component was allowed to have (1) 
an extended width along ?/a which is required to describe 
similar angular correlation data from Au-Au collisions at 
RHIC, (2) non-Gaussian distortions which are suggested 
by possible higher-order fragmentation processes [2j]> [22j , 
and (3) additional soft radiated gluon - jet particle corre- 
lations [3l|, HH . None of these extensions enabled the jet 
related components to describe the same-side correlation 
ridge. 

However, the same-side ridge was well described when 
either an additional same-side Gaussian ridge on az- 
imuth or an azimuth quadrupole [cos(20a)] was in- 
cluded in the fitting. The former structure was predicted 
for dense, strongly interacting systems where initial en- 
ergy/momentum density fluctuations are propagated to 
the final hadronic state via pressure driven radial flow. 
The latter (quadrupole) has been observed in angular 
correlations for relativistic heavy ion collisions. 

Based on the present analysis neither the quadrupole 
or same-side Gaussian ridge phenomenological descrip- 
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tions can be ruled out, although the quadrupole model 
results in smaller % 2 /dof and residuals for both the lower 
p t {p t > 0.1 GeV/c) and higher p t (1 < p t < 3 GeV/c) 
high event multiplicity data. No evidence for either the 
quadrupole or same-side Gaussian ridge was found in 
this analysis for the 7 TeV minimum-bias p-p correla- 
tion data. Alternative mechanisms for producing the 
quadrupole correlation have been proposed [9l-lll| which 
do not rely on rapid thermal equilibration and pressure 
driven collective flow. It would be interesting to apply 
those models to the CMS p-p correlation data. 

The present analysis shows that the CMS ridge could 
be a manifestation of the quadrupole correlation, well- 
known from relativistic heavy ion collisions, and possibly 
appearing in p-p collisions for the first time in these data 
from CMS. The phenomenological results presented here 



and in [26|, 127J together with the problematic assumption 
of rapid equilibration in p-p collision systems, warrant 
further study of the underlying mechanism(s) producing 
the quadrupole correlation in high energy collisions. 
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